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We present precision measurements of the 3 He + 4 He — » 7 Be reaction in the range E c . m . = 0.33 
to 1.23 MeV using a small gas cell and detection of both prompt 7 rays and 7 Be activity. Our 
prompt and activity measurements are in good agreement within the experimental uncertainty of 
several percent. We find •S'(O) = 0.595 ±0.018 keVb from fits of the Kajino theory to our data. We 
compare our results with published measurements, and we discuss the consequences for Big Bang 
Nucleosynthesis and for solar neutrino flux calculations. 
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I. INTRODUCTION 

The 3 He + 4 He — > 7 Be fusion reaction is an important 
step in the solar p-p chain responsible for producing solar 
neutrinos from 7 Be and 8 B decay. Bahcall and Pinson- 
neault [1] identified it as the solar fusion reaction rate 
most in need of further study, due to its relatively large 
uncertainty. It is also reponsible for essentially all of the 
7 Li produced in the Big Bang, a nuclide whose apparent 
primordial abundance remains a puzzle [2] . 

The 3 He + 4 Hc — » 7 Bc reaction rate has been deter- 
mined previously by detecting the prompt 7 rays [3- 
9] and by counting the 7 Be activity [7, 10, 11] (see 
Fig. 1). The accepted value for the zero-energy 5- factor 
for this reaction is based on a 1998 recommendation of 
5(0) = 0.53 ± 0.05 keVb [12], where the relatively large 
uncertainty stems from an apparent disagreement be- 
tween the results from the two methods. Since then, new 
activity meaurements have been published by a Weiz- 
mann Institute group [13], and the LUNA collaboration 
has presented new activity and prompt results [14, 15]. 

Our experiment was designed to measure both the 
prompt 7 rays and the 7 Be activity produced in the same 
irradiation. The prompt 7 yield was measured with a Ge 
detector at 90° with respect to the beam axis. In order 
to contain the 7 Be activity with a high efficiency and ac- 
curately define the active volume, we used a 3 He gas cell 
target 29.7-mm long with a thin, l-/zm nickel entrance 
foil. Gas pressures of 100 and 200 torr were used. The 
copper beam stopper at the end of the gas cell was re- 
moved after He irradiation and counted together with 
the thin Ta gas cell liner to determine the 7 Be activity. 
We measured the cross section in the range E c . m . = 0.33 
to 1.23 MeV with statistical uncertainties of 3% or less for 
each data point, and an overall systematic uncertainty of 
3% for the activity data and 3.5% for the prompt data. 
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FIG. 1: Energy level diagram. 



II. EXPERIMENTAL PROCEDURE 

A. Experimental apparatus and procedure 

We measured the 3 He + 4 He fusion cross sections 
using the Model FN Tandem Van de Graaff acceler- 
ator at CENPA, the Center for Experimental Nuclear 
Physics and Astrophysics at the University of Washing- 
ton. 4 He + beams of 1.5 — 3.5 MeV were produced in 
an rf-discharge ion source mounted in the terminal of 
the accelerator. The beam current on the target cell 
was typically 450 — 500 nA, which, together with beam 
rastering (see below), was designed to achieve good en- 
trance foil lifetime. Between the switching magnet and 
the 3 He gas target, the beam passed through two mag- 
netic quadrupole doublets, a set of X-Y magnetic ras- 
tering coils and two LN 2 traps. A 7-mm water-cooled 
aperture was located just downstream of the second LN 2 
trap. A second water-cooled 7-mm aperture was located 
17 cm further downstream on a moveable aperture plate, 
shown on the left side of Fig. 2. The beam then passed 
through two 8-mm diameter cleanup collimators, a cylin- 
drical electron suppressor, the l-/^m nickel entrance foil 
mounted on a nickel foil holder (with a 10-mm diameter 
aperture) and into the cell as shown in Fig. 2. The aper- 
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FIG. 2: Gas cell and beam collimation geometry. 



tures and collimators were made from Oxygen-Free High 
Conductivity (OFHC) copper and the beamline and the 
region immediately upstream of the gas cell were pumped 
by cryopumps. 

The beam was tuned by focussing through a 1-mm 
diameter aperture located on the moveable aperture plate 
and backed by a tuning beam stop (see Fig. 2). Then the 
beam was rastered to a square distribution somewhat 
larger than 7x7 mm using 19-Hz and 43-Hz X and Y 
raster frequencies, and sent through the 7-mm diameter 
moveable aperture into the gas cell. 

The cell consisted of an aluminum cylinder with a 
32-mm inside diameter, a 29. 7-mm active length, and 
a 30-mm diameter, 0.25-mm thick OFHC copper beam 
stop soldered onto a copper back plate. The choice of 
OFHC copper for the beam stop was motivated by its 
low prompt background, low Li and B contamination, 
and small backscattering probability for low energy 7 Be 
ions. The back plate was air-cooled, and the cell wall was 
lined with 0.025-mm thick Ta foil in order to catch 7 Be 
recoils that did not end up in the stopper. 

Prompt 7 rays were detected in an N-type Ge detec- 
tor with 100% relative efficiency located at 90° with its 
cryostat front face « 5 cm from the center of the cell. 
A 3.2-mm thick Plexiglass absorber was placed between 
the cell and the detector to ensure that (3 particles in co- 
incidence with 7 rays from calibration sources did not 
penetrate into the Ge. The Ge detector was uncolli- 
mated, and shielded on the sides by 20 cm of Pb and 
on the rear by a partial 10-cm Pb shield. We determined 
the total prompt 3 He( 4 He,7) 7 Be cross section from the 
7o yield and the relatively sharp, isotropic 429-keV sec- 
ondary following 7i emission. We used the 70 photopeak 
centroid to determine the mean reaction energy. 4 He + 
4 He background measurements were made at each bom- 
barding energy in order to determine the prompt 7-ray 
background and provide a check on possible contaminant 
7 Be production. 

After irradiation, the stopper was removed from the 



back plate using a knife, and the delayed 478-keV 7 rays 
from 7 Bc decay in the stopper and the Ta liner were 
counted together in a close geometry using a second (off- 
line) 100% Ge detector in an enclosed 20-cm thick Pb 
shield. 



B. Ge calibration 

The Ge detectors were calibrated using 54 Mn, 60 Co, 
133 Ba, 137 Cs, 88 Y, 113 Sn and 203 Hg sources with activi- 
ties in the range 1 — 7 kBq specified to an accuracy of 
±(0.7 - 1.2)% (1<t) [16]. We also used 7 Be and 24 Na 
sources made by the 10 B(p,a) 7 Be and 23 Na(d,p) 24 Na re- 
actions. Penelope [17] Monte Carlo simulations were 
used to calculate finite source size correction factors for 
the online detector. The source line intensities were de- 
termined from asymmetric Gaussian function fits to the 
photopeaks. The energy calibration of the online detec- 
tor was determined using in situ naturally occuring 7- 
ray lines, particularly the 1460.85- and 2614.55-keV lines 
from 40 K and 208 T1 decay. 

1. Offline (activity) Ge detector 

The offline Ge detector efficiency was measured in two 
ways. First, the 7 Be source strength was determined us- 
ing the online Ge detector together with the calibrated 
sources listed above, at a source-to-detector distance of 
25 cm, large enough that 7 — 7 summing was unimpor- 
tant. Then the 7 Be source was used to calibrate the off- 
line Ge in the close geometry, approximately 2.7 mm from 
the front face of the carbon fiber cryostat window. Sec- 
ond, the 54 Mn, 113 Sn and 203 Hg calibration sources were 
used to calibrate the close geometry. The results of these 
two methods agreed well, with an absolute efficiency de- 
termination at 478 keV of typically 0.1246 ± 0.0015. A 
measured ±0.25-mm longitudinal variation in source and 
stopper placement resulted in an additional ±1% system- 
atic error for the different measurements based on a mea- 
sured efficiency sensitivity of 4% per mm displacement. 
A TRIM [19] calculation indicated that the spatial profile 
of the 7 Be implanted in the stopper was approximately a 
9-mm diameter disk. We included a 0.997 correction fac- 
tor for this effect based on off-axis 7 Be source efficiency 
measurements . 

trim calculations for our beam energies, foil and gas 
thicknesses indicated that > 99% of the 7 Be atoms pro- 
duced by 3 He + 4 Hc — > 7 Be should be implanted in the 
beam stop. Since Be adheres readily to surfaces, any 
remaining 7 Be should have deposited on the interior sur- 
faces of the cell, of which 80% was represented by the 
beam stop plus the Ta liner. 

We counted the 7 Be activity with the stopper mounted 
as described above, and the Ta liner cut up and mounted 
directly behind the stopper. In this geometry the 7 rays 
from the liner were counted with an efficiency reduced 
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FIG. 3: (Color online) Absolute online Ge photopeak detec- 
tion efficiency vs. E-,. Solid points - measurements. Error 
bars are statistical. Curve - fit to the data of the function 
e(-E 7 ) = aE^(l + cE^), with a,b and c parameters. Bottom 
panel: fractional fit residuals. 



by 5% due to the greater source-to-detector distance as 
well as absorption in the stopper and Ta. Separate mea- 
surements for several bombarding energies in which only 
the Ta liner was counted indicated < 2% for the frac- 
tion of 7 Be activity on the liner; thus the reduced effi- 
ciency for counting the liner resulted in a negligible addi- 
tional uncertainty. Our overall systematic uncertainty in 
the efficiency for counting the 7 Be activity was ±1.6%. 
We assumed a 7 Be half-life of (53.22 ± 0.06) d and a 
(10.44 ±0.04)% branching ratio for decay to the 478-keV 
daughter state [21]. 



2. Online (prompt) Ge detector 

Online Ge detector efficiency measurements in the 
range E 1 — 280 to 2754 keV were made using the 
24 Na, 54 Mn, 60 Co, 137 Cs, 88 Y, 113 Sn and 203 Hg sources 
mounted in the center of the gas cell. For the two-line 
sources, summing corrections of up to (5 ± 1)% were 
made based on measured peak/total ratios and computed 
7 — 7 angular correlation coefficients. The results, shown 
in Fig. 3, determine the online detection efficiency to 
±1.6%. The position of the online detector was not ex- 
actly the same for the measurements at different bom- 
barding energies. We assigned an additional ±1% vari- 
able systematic error to its efficiency based on numerous 
repeat efficiency measurements. 

We made additional efficiency measurements with the 
sources mounted inside the cell and displaced along its 
axis, as shown in Fig. 4. Increased absorption causes 




FIG. 4: (Color online) Solid (blue) circles: measured online 
Ge detector efficiency vs. position for a 137 Cs source on the 
axis of the gas cell. Solid (red) squares: Penelope simula- 
tions. The curve is to guide the eye. The center of the cell 
is at 2.63 cm, the simulated end points are located at posi- 
tions near the ends of the cell, and the position scale has an 
arbitrary offset. 



the efficiency to drop near the foil holder and near the 
back plate. The measured dependence of efficiency on 
displacement is reproduced well by Penelope [17] Monte 
Carlo efficiency simulations described below. 

The 70 yield was obtained by summing the counts in 
a background-subtracted window around the full-energy 
peak, and the 429-keV yield was determined by fitting 
the peak shape with the same function used for fitting 
the source lines and the 478-keV line. 



3. Penelope Monte Carlo efficiency calculations 

Prompt 7-rays were emitted from an extended cylin- 
drical volume inside the gas cell. The prompt 70 energy 
distribution was broadened by beam energy loss in the 
gas and by the Doppler shift, and as a result the broad 70 
peak contained a 2 — 3% Compton scattering component. 
In order to account accurately for these effects we carried 
out detailed efficiency simulations using the Penelope 
computer code [17]. 

Fig. 5 shows the detector and gas cell geometry used 
in the Penelope calculations. Since we needed sim- 
ulated 7-ray spectra only for energies near the photo- 
peak, we included only materials that resulted in line-of- 
sight absorption and/or small angle forward scattering 
between the cell and the detector. Detector properties 
were taken from manufacturer's specifications [18]. Ef- 
ficiency simulations for the calibration source energies 
were typically (5 — 7)% higher than our measurements 
for E 7 = 0.4 - 3.0 MeV. 

We modeled the 70 detection efficiency using TRIM to 
calculate the energy and position dependence of 7 rays 
emitted from the gas cell for each bombarding energy. 



4 



Ge crystal 



carbon fiber 
window v 

Plexiglass 





gas cell 



3 cm 



FIG. 5: (Color online) Gas cell and detector geometry used 
in Penelope simulations. The beam enters from the left, 
passing through the collimator mounting tube and electron 
suppressor. The Ge detector is shown with its carbon fiber 
window and Plexiglass absorber. 
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FIG. 6: (Color online) E c . m . = 0.43 MeV (E a = 
1.7 MeV) prompt background-subtracted 7-spectrum (points) 
and Penelope simulation (curve) in the region of the 70 peak. 



The measured Ni foil thickness, 3 He gas pressure and ac- 
celerator beam energy were inputs for these calculations. 
The variation of cross section with energy was included 
by assuming a constant iS-factor over the range of energy 
loss in the cell. The result of these calculations served as 
input for the Penelope simulation. In general the sim- 
ulated lineshapes agreed very well with the data, as can 
be seen for the -E c .m. = 0.43 MeV case shown in Fig. 6. 

In order to minimize our sensitivity to calculations, 
we used Penelope efficiency simulations only for ratios 
that are close to unity. We determined the 70 and 7429 
detection efficiencies e(E^) from the relation 



e(£ 7 ) = e p (E 7 ) 



(1) 



where e p (E 1 ) is the point-source efficiency measured 
at the center of the cell, e' d (E y ) is the Penelope 
distributed-source efficiency, e' p (Ej) is the Penelope 
point-source efficiency and E~ l is the mean 7-ray energy. 
Small (several keV) inaccuracies in the simulated 70 en- 
ergy distribution were taken into account by adjusting 
the window used to sum the simulated yield by the dif- 
ference in the measured and simulated 70 centroids. The 
calculated Penelope efficiency ratios were all within 4% 
of unity. This, together with the fact that the measured 
and simulated point-source efficiencies agree well as de- 
scribed above, suggests a systematic uncertainty on the 
simulated efficiency ratio which is less than the typical 
simulation statistical uncertainty of « 0.4%. 



C. Beam current integration 

We determined the number of a particles entering the 
gas cell by integrating the electrical current measured on 
the gas cell, including foil holder and cell support tube, 
which acted as a Faraday cup as shown in Fig. 2. Beam 
collimation and alignment ensured that the beam passed 
cleanly through the foil holder aperture. Secondary elec- 
tron loss from the foil was suppressed with —600 V ap- 
plied to the electron suppressor. The (negative) current 
on the last (8 mm) collimator just upstream of the sup- 
pressor was always less than 1% in magnitude compared 
to the (positive) cell current. The current integrator was 
checked to a precision of better than 1% with two differ- 
ent precision DC current sources and with a square wave 
generator. 

We also made measurements without foil or gas, in 
which we compared the integrated beam current to the 
integrated beam power determined with a calorimeter 
suitable for relative measurements. We used both pro- 
ton and 4 He + beams at 1.00 MeV, and we compared a 
geometry close to that used for the cross section mea- 
surements with a different Faraday cup geometry. In all 
these comparisons, the ratio of the number of beam par- 
ticles determined from the integrated beam current and 
from the integrated beam power was the same within the 
measurement accuracy of ±1%. 

A potential concern is the possiblity of He + charge- 
changing collisions in the residual gas of the beamlinc 
upstream of the target. Pickup is negligible compared 
to stripping [22] , and the strong defocussing effect of the 
beam transport, which would transmit less than 1% of 
a He ++ beam when tuned for He + ions, indicates that 
only those He ++ ions produced in the 2-m beamline sec- 
tion downstream of the last quadrupole lens may reach 
the target cell. In this region the beamline pressure was 
< 2 x 10~ 6 torr. The stripping cross section on O2 and 
N2 shows a broad maximum at E a 1.0 MeV and a 
peak cross section of « 1.6 x 1CP 16 cm 2 [22]. These pa- 
rameters imply the fraction of He ++ ions reaching the 
target was <0.2%. The rastering would reduce the He ++ 
fraction by an additional factor of «4; hence He + charge- 
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changing collisions were not an important concern for our 
measurements. 



D. Target cell gas and foil properties 

The nickel entrance foils were 1-^tm nominal thickness, 
glued onto the nickel foil holders. Fresh foils were used 
for most bombarding energies. Foil thicknesses were mea- 
sured to ±1% using the 3.2-MeV a particles from 148 Gd. 
Most foils leaked He at a slow rate, and the cell pressure 
was "topped off' in order to maintain a constant average 
pressure over time. Pressure variations over the course 
of a 1-hour irradiation were < 6%. The increase in cell 
depth due to foil bulge was measured to be 0.3 — 0.4 mm 
averaged over the 7 — 8-mm beam diameter. 

We used 3 He gas with 99.999% (99.99%) chemical 
(isotopic) purity, and naturally occuring He gas with 
99.9999% chemical purity. The cell was flushed with 3 Hc 
(or 4 He) before and during the cross section measure- 
ments, which were carried out at nominal pressures of 
200 torr (two highest energies) and 100 torr. Gas purity 
was measured in situ after some time in the cell, using 
a residual gas analyzer and found to > 99% pure. Gas 
pressure was measured with a Baratron [20], whose pre- 
cision was checked with a mercury manometer to better 
than 1% at 100 — 200 torr and was recorded in 1-s time in- 
tervals with an ADC (Analog-to-Digital Converter). Cell 
temperature was measured with an insulated probe at- 
tached to the cell wall and recorded at the beginning and 
end of each one-hour run. 



E. Contaminant r Be production 

As Hilgemeier et al. [9] pointed out, 7 Be may be 
produced by a background process, e.g. 6 Li(c?,n) 7 Be or 
10 B(p,a) 7 Be if there is both a (proton or deuteron) con- 
tamination in the 4 He beam and a ( 6 Li or 10 B) contam- 
inant in the gas cell, e.g. the foil or beam stop. This 
concern was apparently overlooked in earlier 3 He + 4 He 
— > 7 Be activity measurements. We guarded against this 
possibility by counting the stoppers from the 4 He + 4 He 
irradiations for 5 of the 8 bombarding energies. We ob- 
tained 1 — 2% upper limits (la) on contaminant 7 Be pro- 
duction from each of these measurements. 

In order to understand the possible magnitude of this 
problem in older 3 He + 4 Hc — > 7 Bc activity measure- 
ments, we made a study of beam and target stopper 
impurities. Using our 90° analyzing magnet and a 3- 
MeV 4 He + beam, we found a and/or DH^ satellite 
beam with intensity 0.1 — 1% (depending on source prepa- 
ration) at slightly higher rigidity than the 4 He + beam. 
Accelerator voltage instabilities, e.g. sparks, could allow 
this beam to pass through the analyzing magnet slits 
with a non-zero duty factor. We then bombarded vari- 
ous Co, Ni, Cu, Nb, Ta, Pt and Au materials selected for 
good chemical purity (> 99.9%) with 0.75-MeV protons 



TABLE I: Beam heating density correction factors, percent 
per 600 nA. 



Case 


£ rcs (McV) 


100 torr 


200 torr 


M Mg(a, 7 r 8 Si 


3.20 




3.9±0.6 


10 B(a,p- 7 ) 13 C 


1.51 


5.3±0.8 


7.8±1.0 



and 1.5-MeV deuterons (the energies that would result in 
the above example) and measured the 7 Be activity pro- 
duced [23]. Our results suggest that contaminant 7 Be 
production in the older activity experiments was proba- 
bly not significant. We note that the Osborne et al. activ- 
ity measurements [7] were made with a Tandem acceler- 
ator in which the beam was stripped in the terminal and 
hence was unlikely to contain molecular impurities. We 
also note that a 3 He beam experiment is more danger- 
ous since 7 Li( 3 He,t) and 6 Li( 3 He,rf) reactions can lead 
to contaminant Be production without the need for a 
beam contaminant. 



F. Beam heating corrections 

The region of the gas cell illuminated by the beam has 
a higher temperature and hence lower density than the 
rest of the cell. We measured this effect using narrow 
resonances in the 24 Mg(a, 7) 28 Si and 10 B(a,p - 7) 13 C 
reactions at E [cs — 3.20 MeV and 1.51 MeV, respectively. 
We coated the gas-cell stopper with a thin layer of natural 
Mg or B and measured the thick-target resonance yield 
with and without the foil and gas. Then, with the foil 
and gas present, and the beam energy adjusted halfway 
up the leading edge of the resonance, we measured the 
resonance yield as a function of the beam current. The 
gas density correction factors determined from these data 
are shown in Table I. The beam heating corrections for 
our 200-torr 3 Hc + 4 Hc — > 7 Be measurements were de- 
termined by interpolating between the entries shown in 
Table I, while the 100-torr 3 He + 4 He — > 7 Be corrections 
were determined by extrapolating the 100-torr measure- 
ment, scaling by the calculated gas dE/dx— see Table II. 

G. 70 anisotropy 

The 429-kcV 7 ray following 71 emission is isotropic, 
since it proceeds from a J — 1/2 state; however, 70 need 
not be. Tombrello and Parker [24] calculated the angular 
anisotropy, which they found to be small but nonzero in 
our energy range. Column 6 of Table III shows the 90° 
70 anisotropy correction factors 1/(1 — a 2 Q2/2) calcu- 
lated from [24], where a 2 is the coefficient of the Legen- 
dre Polynomial P 2 and Q 2 = 0.9 is the angular attenua- 
tion coefficient for our geometry. These corrections range 
from 0.99 to 1.04 for our energies. Kim et al. [25] calcu- 
lated the 70 angular distribution at several energies; at 
^c.m. =0.1 MeV their result agrees well with [24] while at 
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FIG. 7: (Color online) Prompt 7-spectra measured at E a — 
1.7 MeV (E c .m. = 0.43 MeV). Top panel, from top to bot- 
tom: 4 He + 3 He (black), 4 He + 4 He (red), and beam-off 
background (blue) spectra. Bottom panels: background sub- 
tracted spectra in the region of the 429 keV line and 70 . The 
window used to determine the 70 yield is indicated. 



E c .m. — 1 MeV their result corresponds to an anisotropy 
correction factor 2% larger than [24]. 

The angular distribution is sensitive to the phase dif- 
ference between s- and d-wave capture amplitudes, which 
has not been measured in our energy range. Both calcula- 
tions assumed a hard-sphere extrapolation of the nuclear 
part from data at higher energies, where measurements 
are consistent at the level of 5° [26]. We estimate that 
a 5° change in this phase difference would give rise to a 
2% change in the anisotropy correction at our energies, 
on average. Thus it seems reasonable to assume a ± 2% 
uncertainty in this correction. We note that Krawinkel et 
al. [6] are the only ones to have measured the anisotropy, 
and their results do not constrain it at this level of pre- 
cision. 



III. RESULTS 

We measured the 3 He + 4 He — > 7 Bc fusion cross section 
at eight bombarding energies from E c _ m . — 0.33 to 1.23 
MeV. Some properties of these measurements are shown 
in Table II. At each bombarding energy we measured the 
4 He + 3 He yield in a series of 1-hour runs for the total 
time indicated. The rastered beam current on the cell 
was 420 - 500 nA. 
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FIG. 8: (Color online) Prompt 7-spectra measured at E a = 
3.5 MeV (£ c .m. = 1.23 MeV). The colors correspond to those 
used in Fig. 7. See Fig. 7 caption for additional information. 



4 He + 4 He background measurements were also made 
at each bombarding energy, for an irradiation time 60 — 
130% of the 4 He + 3 He irradiation time. Some back- 
ground increase with time was observed for most of the 
runs, which we attributed mainly to the 13 C(a,n) 16 re- 
action on carbon buildup. 4 He + He measurements were 
made before and after the 4 He + 3 He measurements, and 
summed, in order to correct, to first order, for the prompt 
7-ray background increase with irradiation time. Sepa- 
rate stoppers were used for the 4 He -I- 3 He and 4 Hc + 4 Hc 
irradiations so we could count the 4 He + 4 He stoppers to 
check for contaminant 7 Be production (see Sec. HE). 

The 7-ray spectra measured at E c . m . = 1.23 and 0.43 
MeV are shown in Figs. 7 and 8. For most of the mea- 
surements the beam-related background (per /iC) deter- 
mined from 4 He + 4 He was statistically consistent with 
the background observed above the 70 peak in the 4 He 
+ 3 He spectra. Even so, we adopted a more conservative 
procedure of determining the 4 He + 3 He background by 
normalizing the (beam-related part of the) 4 He + 4 Hc 
spectrum to a region above the 70 peak in the 4 He + 
3 He spectrum. The 4 He + 4 He spectra shown in the 
top panels of these figures are shown with this normal- 
ization. The lower panels in Figs. 7 and 8 show the 
background-subtracted prompt spectra in the region of 
the 70 and 429-keV peaks, and the window used to sum 
the 70 counts. Sharp features in the subtracted spec- 
tra near 2614 and 1765 keV originate from large back- 
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TABLE II: Conditions of the cross section measurements. 



£ c . m .(MeV) 


E a (MeV) 


Foil(/im) 


A£ foil (MeV) 


P gas (torr)° 


A£ gas (MeV) 


Iav(nA) 


BH(%) 6 


T irrad (hr) c 


T act (d) 


0.3274 


1.5 


0.968 


0.677 


100 


0.166 


480 


5.3 


102.3 


26.6 


0.4260 


1.7 


0.937 


0.641 


100 


0.158 


460 


4.9 


116.5 


12.0 


0.5180 


1.9 


0.955 


0.638 


100 


0.148 


450 


4.5 


30.9 


33.6 


0.5815 


2.1 


1.069 


0.691 


100 


0.141 


420 


3.9 


24.4 


20.5 


0.7024 


2.35 


1.069 


0.665 


100 


0.129 


450 


3.9 


23.3 


10.4 


0.7968 


2.6 


1.072 


0.642 


200 


0.240 


500 


5.9 


20.3 


7.2 


1.2337 


3.5 


1.065 


0.545 


200 


0.182 


440 


3.5 


14.7 


6.8 


1.2347 


3.5 


1.051 


0.538 


200 


0.182 


120/500 


1.2/4.9 


41.8 


31.8 



"Nominal gas pressure. 
'Average beam heating correction. 
c4 He + 3 He irradiation time. 
^Activity counting time. 
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FIG. 9: (Color online) E c . m . = 1.23 and 0.43 MeV {E a = 
3.5 and 1.7 MeV) offline activity spectra. Top: 4 He + 3 He 
irradiation. Middle: 4 He + 4 He irradiation. Bottom: source- 
out background. Counting periods (in d) are shown. 



ground peaks that have been subtracted imperfectly (due 
to small gain shifts) and have small net areas. 

In our i? c . m . = 0.80 MeV data we observed a large 
carbon buildup. This allowed us to check that the beam- 
related backgrounds due to carbon buildup and the back- 
ground present before buildup had the same shape in the 
region of interest. Thus our method of background deter- 
mination under the prompt 70 peak appears to be robust. 

The E CAn ,— 1.23 and 0.43 MeV activity spectra are 
shown in Fig. 9. Neither the source-out background 
(measured with a Cu stopper that had not been irradi- 
ated) nor the 4 He + 4 He stopper measurements showed 
any sign of a 478-keV peak. 

We determined the mean reaction energy from the 70 
centroid. The advantage of this method is that it does 
not depend on knowledge of the beam energy or the en- 
ergy loss in the foil or gas, and it is determined directly 
from the events used to determine the cross section. How- 



ever, the broad 70 peak contains a small forward-angle 
Compton scattering component in addition to the major- 
ity of events in which the full 7-ray energy is deposited 
in the Ge. The 70 centroid must be corrected for this 
component in order to obtain the mean reaction energy. 
(Note that the Compton contribution to the 70 detection 
efficiency is included in the Penelope ratio factor that 
appears in Eq. 1). We computed centroid corrections of 
1.2 ± 0.1 and 1.9 ± 0.1 keV for the E c . m = 0.33 and 
1.23 MeV cases using the measured detector response for 
£J 7 = 1836 and 2754 keV, and we interpolated linearly 
to obtain the corrections for the other reaction energies. 
These corrections shift our inferred ^-factors down by 
1.2% (0.9%) for the E c , m = 0.33 (0.43) MeV cases and 
less at the higher energies. 



We also corrected for a small net positive Doppler shift 
in the 70 centroid. This resulted primarily because, on 
average, the 7-rays were emitted somewhat (< 2 mm) 
upstream of the center of the cell. The calculated cor- 
rection is 0.35 (0.2) keV for the E c . m . = 0.33(0.43) MeV 
cases, and has the opposite sign to the Compton correc- 
tion discussed above. All of our quoted -E c .m. values and 
5-factors include these corrections. 



As an additional check on beam current integration, 
we remeasured the cross section at E CAn , — 1.23 MeV us- 
ing 4-mm instead of 7-mm diameter beam-defining aper- 
tures. This further minimized the secondary electron 
currents on the cleanup collimators. We determined the 
beam heating correction for this collimation geometry 
directly from 3 He + 4 He — > 7 Be data measured at 120 
and 500 nA. The good agreement between the activity S- 
factors deduced from these two E c , m , = 1.23 MeV cross 
section measurements provides additional evidence that 
our beam current integration and beam heating estimates 
are reliable. The prompt ^-factors for the two 1.23 MeV 
measurements agree less well; however, here we measured 
the 4 He -I- 4 He beam-related background only at 500 nA. 
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TABLE III: Our ^-factors, prompt 7-ray branching ratios, and calculated 70 anisotropy. 



£c. m .(MeV) & ct (keVb) a 5prom.pt (keVb)"-'' S prompt /Szct b 7429/70 b Anisotropy^ 

0.3274 ± 0.0013 0.495 ± 0.015 0.492 ± 0.014 0.994 ± 0.042 0.410 ± 0.023 0.988 

0.4260 ± 0.0004 0.458 ± 0.010 0.438 ± 0.006 0.956 ± 0.024 0.405 ± 0.009 1.002 

0.5180 ± 0.0005 0.440 ± 0.010 0.421 ± 0.007 0.957 ± 0.026 0.394 ± 0.009 1.013 

0.5815 ± 0.0008 0.400 ± 0.011 0.398 ± 0.008 0.995 ± 0.035 0.422 ± 0.013 1.019 

0.7024 ± 0.0006 0.375 ± 0.010 0.382 ± 0.006 1.020 ± 0.031 0.424 ± 0.010 1.030 

0.7968 ± 0.0003 0.363 ± 0.007 0.371 ± 0.004 1.022 ± 0.023 0.427 ± 0.005 1.036 

1.2337 ± 0.0003 0.330 ± 0.006 0.327 ± 0.004 0.990 ± 0.021 0.439 ± 0.006 1.040 

1.2347 ± 0.0003 0.324 ± 0.006 0.340 ± 0.004 1.050 ± 0.023 0.443 ± 0.007 1.040 



"Uncertainties include statistics, ±1% variable systematic and the 
contribution from the -B c .m. uncertainty. 
''Assuming 70 isotropic. 

c 70 anisotropy correction 1/(1 — 112Q2/2) calculated from [24]. 



TABLE IV: Systematic uncertainties AS , (£ c . m .)/S'(S c . m .). 
The variable uncertainties are assumed random between the 
various S(E c . m .) values, while the scale factor uncertainties 
are the same for all the values. 





Type" 


AS/S(%) 


Origin 


eV b 


variable: 


a, p 


1.0 


source positioning 




scale factor: 


a, p 


1.6 


Ge efficiency b 


E 




P 


1.0 


beam centering 


6 
LU 




P 


1.5 


70 anisotropy 


Tj 




a, p 


1.6 


current integration 


c/o 




a, p 


1.3 


gas pressure 






a, p 


0.6 


temperature 






a, p 


0.7 


gas purity 






a, p 


0.4 


gas cell length c 






a, p 


1.0 


beam heating 




total scale factor: 


a 


3.0 




"to 




P 


3.5 




"O 

ifl 




common 


2.7 




CD 



"Activity and/or prompt. 

'Prompt value includes summing and simulation uncertainties. 
Includes foil bulge. 




0.04 


-0.04 









t i 



0.2 0.4 



0.6 
E 



0.8 
c m. (MeV) 



1.2 1.4 



.S'-factors 



FIG. 10: (Color online) Solid (blue) points: .Sact^c.m.) vs. 
-Bc.m.. Curve: fit of Kajino's function. Bottom panel: frac- 
tional fit residuals. 



Our measured S'-factors are shown in Table III. The 
mean center-of-mass reaction energy was determined by 
subtracting the 1586.1 ±0.1 keV Q-value [27] from the 
70 centroid. Cross sections were converted to S'-factors 
using the relation 



S{E c . m .) — <r(E c . m .)E c 



,(E G /Bo. m .) 1/2 



(2) 



with Eq /2 = 164.13 keV 1 / 2 . The tabulated uncertain- 
ties in the S'-factors and their ratios include statistics, 
the ±1% variable systematic uncertainties shown in Ta- 
ble IV, and the contribution of the uncertainty in the 
mean reaction energy E c m . - the latter is important only 
at the lowest bombarding energy. We also show the ex- 
perimental branching ratio 7429/70 ( = 71/70) an d the 
anisotropy calculated from [24] (see Sec. II G). 

Table IV shows the sources of systematic uncertainty 



and their contributions to the S-factor uncertainty. The 
±1% variable systematic uncertainties given in Table IV 
are included in the S-factor values and ratios shown in 
Table III and in the fitting of S(E c . m .) vs. E c . m ., while 
the systematic scale-factor uncertainty is folded with the 
fit uncertainty to give the tota l unce rtainty. The fit un- 
certainties include the factor \J x 2 l v whenever x 2 l v > 1- 

Fig. 10 shows S , ac t(£'c.m.) vs. E c , m , together with a 
fit in which the theoretical shape given by Eq. 6 of Ka- 
jino et al. [28] has been scaled to fit the data. Fig. 11 
shows 5 , p rom pt(-E ; c.m.) based on an assumed isotropic an- 
gular distribution for 70 together with the Kajino fit. In 
these fits the shape parameters of the fit function arc held 
constant and only one parameter, S*(0), is varied. 

Fig. 12 shows the ratios S pmmp t/S ac t for both the 
isotropic and anisotropic assumptions for 70. We fit 
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TABLE V: Our 5(0) values and prompt/activity ratios. 
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FIG. 11: (Color online) Solid (blue) points: 5 pr ompt(-Ec.m.) vs. 
-Bc.m., for the isotropic 70 assumption. Curve: fit of Kajino's 
function. Bottom panel: fractional fit residuals. 
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FIG. 12: (Color online) 5p r0 mpt /S ac t vs. E c . m .. Solid (blue) 
circles: assuming 70 isotropic. (Red) x's: 70 corrected for 
theoretical anisotropy as given in Table III. The crosses have 
been displaced from the circles for clarity. 



the ratios S pr0 mpt /Sact vs. E c , m , with a horizontal line, 
and we also fit our S promp t(-E c .m.) values based on the 
anisotropic assumption. The results of all these fits are 
shown in Table V, including the confidence level P(x 2 > v) 
for each fit. The fit to S ac t(-E c .m.) has a good qual- 
ity. The fit to Sp r ompt (-Ecm.) is very poor quality for 
the anisotropic case, and although the confidence level is 
still low, it's much better for the isotropic assumption. 
The ratio fits are also poorer for the anisotropic assump- 
tion, and acceptable for the isotropic assumption. Thus 
we conclude that the 70 anisotropies calculated from [24] 
are not consistent with our data, while the isotropic ap- 
proximation is better. Thus our best S(0) values based 



Quantity 


Value" 


x> 6 


p(%r 


Sact(0) 


0.595 ±0.018 keVb 


0.72 


65 


^prompt (0) 


0.596 ±0.021 keVb 


2.7 


1 


^prompt /Sact 


0.999 ± 0.030 


1.8 


9 


'S'prompt (0) 


0.607 ±0.022 keVb 


5.3 


0.0004 


'S'prompt / Sact 


1.015 ±0.032 


2.9 


0.5 


Combined 3(0Y 


0.595 ± 0.018 keVb 







"Uncertainties include systematic contributions. 
b u = 7. 

C P(X 2 , v) is the confidence level of the fit. 
d 70 assumed isotropic. 



e 70 corrected by theoretical anisotropy 
f Activity and isotropic prompt. 



see Table III. 



(3) 



on the Kajino fit function are 

Sact(O) = 0.595 ± 0.018 keVb, 
and 

^prompt (0) = 0.596 ± 0.021 keVb. (4) 

Our best value for 5(0) is obtained by combining the 
activity and prompt results given above, taking into ac- 
count the 2.7% common scale factor uncertainty given in 
Table IV: 



IV. 



5(0) = 0.595 ±0.018 keVb. 



COMPARISON WITH OTHER 
EXPERIMENTS 



(5) 



We found that fitting a given data set with different 
theories resulted in 5(0) values that varied typically by 
several percent (5% in the case of ref. [10], which con- 
sists of one high-energy point). Hence in comparisons of 
different experimental results, it is important to use the 
same fitting function - the exceptions are the LUNA re- 
sults [14, 15] which were determined from data at very 
low energy. Since the resonating-group calculations of 
Kajino et al. do a good job of reproducing other ob- 
servables, and result in a reasonable theoretical value 
for 5(0), we chose this fit function. We obtained data 
from the original publications, from the NACRE compi- 
lation [29] and from other sources (see Refs. [4, 7]). 

The results are shown in Table VI. Most differences 
between the earlier results shown here and in Table II 
of [12] are due to the differing fitting functions; the no- 
table exception is the Osborne et al. [7] activity result, 
which is significantly higher here. This problem has been 
noted before - see Kajino et al. [28] Table 1 and footnotes. 
We fit Osborne et al.'s S ac t(E c , m .) values quoted in their 
published Table 1 [7], which are the same as shown in 
their published Fig. 9 and in Osborne's Ph.D. thesis [7]. 
Volk et al. [11] give insufficient detail for us to interpret 
their integral measurement in a manner consistent with 
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TABLE VI: S(O) values from our work and from published 
data. 





S(0) (keVb) 


Reference 




activity: 


0.577 ± 0.035 


Osborne et al. 


[7] 




0.660 ± 0.040 


Robertson et al. 


[10] 




0.560 ± 0.030" 


Volk et al. 


[11] 




0.546 c ± 0.020 


Nara Singh et al. 


[13] 




0.545 ± 0.017 


Gyurky et al. 


[14] 




0.595 ± 0.018 


Present work 




prompt: 


0.481 ± 0.053 


Parker et al. 


[4] 




0.579 ± 0.07 d 


Nagatani et al. 


[5] 




0.449 ± 0.06 E 


Krawinkel et al. 


[6] 




0.522 ± 0.03 


Osborne et al. 


['] 




0.478 ± 0.04 


Alexander et al. 


[8] 




0.542 ± 0.03 


Hilgemeier et al. 


[9] 




0.560 ± 0.021 


Confortola et al. 


[15] 




0.596 ± 0.021 


Present work 




total:-' 


0.560 ± 0.017" 


Confortola et al. 


[15] 




0.595 ± 0.018 


Present work 





"From our fits of Kajino's function to published S(E) values, with 
uncertainties as quoted by authors except where noted. 

"Value and uncertainty as quoted by authors. 

c Scaled by 1.008 to account for the different 7 Be decay BR as- 
sumed. 

^Uncertainty as quoted in [12]. 

e Value corrected by xl.4, and uncertainty as given in [9]. 
^Combined activity and prompt. 
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f Robertson 


^ Present 


I Volk 





— Prompt fit 


J Alexander 
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FIG. 13: (Color online) Left panel, (red) points: S'act(O). 
Right panel, (blue) points, 5 , pr ompt(0) . The horizontal lines 
are the fits to the activity values (left panel) and prompt val- 
ues (right panel) as given in Table VII. 



our other analyses; hence we simply quote their published 
value. 

The activity and prompt 5(0) values from Table VI are 
shown in Fig. 13 together with separate fits to the activ- 
ity values and to the prompt values. These and other fit 
results are shown in Table VII. Our contaminant study 
reported in Sec. HE suggests earlier activity results prob- 
ably did not suffer from contaminant 7 Bc production. 
Earlier theoretical work showed that processes other than 
single-photon emission, such as E0 emission, should be 



— Fit to all 


except Gyurky 


^ Nara Singh 


J Confortola total 


J Gyurky 


J Present total 



0.65 



0.60 - 



0.55 



0.50 



FIG. 14: (Color online) Recently measured S(0) values. 
Points, from left to right: Ref. [13] and Ref. [14] activity val- 
ues, Ref. [15] and present activity+prompt values. Horizontal 
line - fit to Ref. [13], Ref. [15] total and present total as given 
in Table VII. 




TABLE VII: S(0) comparisons. 



Data 



all activity 

all activity except [10] 

all prompt 

all c 

present tot ± [15] tot ± 
present tot + [14] + [13] 



[13" 



S(0) a 
0.568 ±0.014 
0.563 ±0.011 
0.549 ±0.016 
0.559 ±0.012 
0.568 ±0.014 
0.562 ±0.017 



x 2 A 



p(%)" 



2.1 
1.3 
2.0 
2.1 
1.8 
2.5 



5 
4 
7 
11 

2 
2 



6 
27 
5 
2 
16 



"All fit uncertainties contain the factor \Zx 2 / u 

"P(xV)- 

C A11 except the last 2 prompt values and the last 2 activity values 
in Table VI. 



negligible [30], while our measured prompt/activity ra- 
tios imply an experimental uppper limit on E0 emission 
of «2%. Thus we see no reason why older activity mea- 
surements, as a group, should be flawed. However, as 
shown in Table VII, the fit to all activity values has a 
poor x 2 l v an< i a poor confidence level, though an ac- 
ceptable quality fit is obtained if the Robertson et al. 
datum [10] is removed. 

Comparing recent measurements, our total 
(prompt±activity) S(0) agrees with that of Con- 
fortola et al. [15] (see Table VI), and is higher than the 
activity results of refs. [13, 14], as shown in Fig. 14. 
Combining our total value with that of [15] and with the 
activity result of [13] yields a minimally acceptable fit; 
substituting the activity result of [14] for the [15] total 
value yields a similar result though the fit quality is not 
as good - see Table VII. Note that we do not include the 
results of both [14] and [15] in the same fit since these 
values are presumably highly correlated, and the degree 
of correlation is not specified. We also did not derive an 
activity value for S(0) from [15] for the same reason. 
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V. CONCLUSION 

We have presented new results for the 3 He + 4 He — ► 
7 Be 5-factor based on measurements of the prompt 7- 
rays and the Be activity produced in the same irradia- 
tion. We find good agreement between the two methods, 
and a combined result 

5(0) = 0.595 ±0.018 keVb, (6) 

based on fits of the theory of Kajino et al. [28] to our 
data. 

The lack of good agreement between the various 5(0) 
values precludes a meaningful determination of a "best" 
value. However, the recent, precise 5(0) determina- 
tions (ours and refs. [13-15]) all lie between 0.53 and 
0.61 keVb, suggesting that the true value of 5(0) also 
lies in this range. Comparing to older recommenda- 
tions made before the recent measurements, this range 
is is well inside the range recommended by NACRE [29] , 
0.54 ± 0.09 keVb, and somewhat higher than the ranges 
recommended by [12], 0.53 ± 0.05 keVb, and by [31], 
0.51 ± 0.04 keVb. 

Further progress will probably depend on new mea- 
surements. It may be important to note that a new mea- 
surement using a recoil separator to detect the 7 Be nuclei, 
previewed in [32], is in good agreement with the present 
results in the region of overlap near £ c . m . = 1.2 MeV. 

VI. IMPACT ON SOLAR MODEL 
CALCULATIONS AND BIG BANG 
NUCLEOSYNTHESIS 

Our value for 5(0) quoted above is 12% larger than the 
accepted value quoted in [12]. What would be the con- 
sequence of an upward adjustment of this size for solar 



model neutrino flux calculations and Big Bang Nucle- 
osynthesis (BBN)? 

It would increase the calculated Standard Solar Model 
(SSM) neutrino fluxes from decay of 7 Be and 8 B in the 
sun by approximately 10% [1]. Using old (new) heavy ele- 
ment abundances, Bahcall et al. [33] quote 1.09 (0.87) for 
the ratio of the calculated SSM to measured 8 B neutrino 
flux, where the measured flux has ± 9% uncertainty [34] 
and the SSM flux has ± 16% uncertainty (both la). Thus 
a 10% upward adjustment of these calculated 8 B fluxes 
would not lead in either case to a significant disagreement 
with the measured flux, based on these uncertainties. 

The uncertainty in the 3 He + 4 Hc -> 7 Bc 5(0) is one 
of the largest uncertainties that contribute to the SSM 
neutrino flux from 7 Be decay in the sun [1]; hence an 
eventual reduction in this uncertainty will reduce signif- 
icantly the uncertainty in the calculated SSM flux that 
the BOREXINO [35] and KamLAND [36] experiments 
are designed to observe. 

The 7 Li problem in Big Bang Nucleosynthesis is well- 
known and long-standing: 7 Li, which results mainly from 
decay of 7 Bc produced by the 3 He + 4 He — ► 7 Be reac- 
tion, is overpredicted by BBN calculations [2]. Thus an 
increase in the 3 He + 4 He — ► 7 Be reaction rate at BBN 
energies would worsen this discrepancy. 
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